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Abstract 

We study the small oscillations regime (RPA approximation) of the 
time-dependent mean-field equations, obtained in a previous work |IJ, 
which describe the time evolution of one-body dynamical variables of a 
uniform Chiral Gross-Neveu system. In this approximation we obtain 
an analytical solution for the time evolution of the one-body dynam- 
ical variables. The two-fermion physics can be explored through this 
solution. The condition for the existence of bound states is examined. 



PACS Numbers : 03.65.Ge, 03.65.Nk, 11.30Rd, ll.lOKk 

1 Introduction 



In a previous work [Ij we obtained in Gaussian mean-field approximation the 
effective dynamics of one-fermion and pairing densities of a off-equilibrium 
spatially uniform (1+1) dimensional self- interacting fermion system described 
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by Chiral Gross- Neveu model (CGNM) JJ]. These dynamical equations ac- 
quire the structure of (colisionless) kinetic equations. They determine the 
time evolution of the all one-fermion densities of this system for a given ini- 
tial condition. Spatial uniformity (translation and reflection invariance) is 
assumed in our derivation. 

Studing the static solutions of these equations in order to renormalize the 
theory IJ, we found an effective potential similar to that obtained by Gross 
and Neveu using the 1/N expansion 0. We also showed that other static 
results which have been discussed in the literature |2|, [| £§] such as dynam- 
ical mass generation due to chiral symmetry breaking and a phenomenon 
analogous to dimensional transmutation can be retrieved from this formula- 
tion. Finally, in |TJ], we obtained numerical solutions for the time evolution of 
the one-body dynamical variables initially displaced from equilibrium. The 
time evolution of the symmetrical and broken chiral phases of our system are 
discussed. 

In this work we explore a particular application of the renormalized non- 
linear mean- field equations obtained in jjj]. We follow the recent work of 
Kerman and Lin || in order to study the near equilibrium dynamics around 
the stationary solution as a tool to investigate the two-fermion dynamics. 
In particular, the resulting equations can be solved analytically to reveal a 
two-(quasi)fermion bound state solution. 

This paper is organized as follows. In Sec. II we linearize the mean-field dy- 
namical equations which describe the time evolution of a off-equilibrium spa- 
tially uniform (1+1) dimensional self-interacting fermion system described 
by Chiral Gross- Neveu model (CGNM). A self-consistent renormalization 
scheme is necessary [I], ||. In Sec. Ill, making use of an analogy with scatter- 
ing theory [0, we obtain a closed analytical solution for the time evolution 
of the one-fermion densities in this regime. Studyng the two-fermion physics 
in Sec. IV, we find the condition for the existence of bound states. Finally, 
the Sec.V is devoted to a final discussion and conclusions. 
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2 Mean-field kinetic equations 



We begin this section by reviewing our approach which describe a formal 
treatment of the kinetics of a self-interacting quantum field. This approach 
was developed earlier for the nonrelativistic nuclear many-body dynamics by 
Nemes and de Toledo Piza || and was more recently applied in the quantum- 
field theoretical context to the self-interacting A0 4 theory in (1+1) dimensions 
The general ideia is to focus on the time evolution of the one-fermion 
and pairing densities. These observables are kept under direct control when 
one works variationally using a Gaussian functional ansatz and will therefore 
be refered to as Gaussian observables. 

We consider an off-equilibrium, spatially uniform, (1+1) dimensional sys- 
tem of relativistic, self-interacting fermions described by Chiral Gross-Neveu 
model (CGNM) M. The Hamiltonian density is given by 
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where £ is a constant which indicates whether the model is invariant under 
discrete 75 transformation (£ = 0) or under the Abelian chiral U(l) group 
(£ = 1). In the form considered here, this is a massless fermion theory in 
(1+1) dimensions with quartic interaction. The model contains N species 
of fermions coupled symmetrically. In the Heisenberg picture, the if} 1 are 
complex Dirac spinors 
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where b^ l and &k.i 2 an d ^,2] are fermion creation and annihilation opera- 
tors associated to positive [negativej-energy solution Ui(k) [w2(k)] of Dirac's 
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equation. 

This model is essentially equivalent to the Nambu-Jona-Lasinio model 
[|IIJ, except for the fact that in (1+1) dimensions it is renormalizable. More- 
over, it is one of the very few known field theories which are assymptotically 
free. To leading order in 1/N expansion ||, the CGNM exhibits a number 
of interesting phenomena, like spontaneous symmetry breaking, dynamical 
fermion mass generation and dimensional transmutation. 

The state of this system (assumed spatially uniform) is given in terms 
of a many-body density operator T of unit trace. Our implementation of 
the Gaussian mean-field approximation consists in approximating this ob- 
ject by a truncated many-body density operator J-'oi't), also of unit trace, 
written as the most general hermitian Gaussian functional of the field opera- 



tors consistent with the assumed uniformity of the system [11]. It will thus be 



written as the exponential of a general quadratic form in the field operators, 
which can be reduced to diagonal form by suitable canonical transformation. 
The most general transformation would in general break both the chiral and 
charge symmetries of the CGNM. In the following development we restrict 
ourselves for simplicity to a special class of transformations (to be called 
Nambu transformations) which break the chiral symmetry only. They can 
be parametrized in a form that incorporates the unitarity constraints as 
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where reflection symmetry of the uniform system is implemented by making 
the parameters (p^ and 7k dependent only on the magnitude of k. 

The Gaussian truncated density operator JF (t) acquires a particularly 
simple form when expressed in terms of the Nambu quasi-fermion operators 
which diagonalize the associated quadratic form, namely 
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Mt) = II [^Ut)f3 K x(t) + (1 - ^,x)Mt)f3l x (t)} , (4) 



k.A 



where z/^a for A = 1, 2 are the Nambu (quasi-fermion) occupation numbers. 

With the help of Eq.(Q) it is an easy task to express ^(x) and ip(x), Eq.(||), 
in terms of /3 kA (t) and Pk,x(t) for A = 1,2. In doing so, one finds that the 
plane waves of ip(x) and if>(x) are modified by a complex, moment- dependent 
redefinition of m involving the Nambu parameters y?k(£) and 7k(t)- The 
complex character of these parameters is actually crucial in dynamical situ- 
ations, where the imaginary parts will allow for the description of time-odd 
(velocity-like) properties. Finally, the mean values of the Gaussian observ- 
ables are parametrized in terms of the (pk{t) and 7k(£) and of the occupation 

numbers z/ k , A (t) = Tr [/^(O/M*)-^*)] for A = 1,2. 

The next step is to obtain the mean-field time evolution for the mean 
values of the Gaussian observables in the context of the initial-value problem. 
In other words, we want the Gaussian mean-field equations of motion for 
the Nambu parameters (pk(t) , 7k(£) and for the quasi-particle occupation 
numbers v^xif). In Ref.jlj] we obtained 

z> k)1 = and z> kj2 = ; (5) 

r. . . . . i _ i7 . Tr ([/3-k,l/3k,2, ^CGNmJ^o) ra s 

[up k + 7 k sin v?k cos y?kj e 7k = r . (6 ) 

Eq.(U) shows that the occupation numbers of the Nambu quasi-particles 
are constant. This is the general isoentropic character of the mean-field 



approximation JT2|, |13 |. The complex equation of motion (y) describes the 
time evolution of the Nambu parameters. From the right-hand side of the 
Eq.(|6|), we see that to obtain the time evolution of the Nambu parameters, 
we have to express the CGNM Hamiltonian in the Nambu basis. 

From the Hamiltonian density ([J) we can explicitly evaluate the Hamil- 
tonian of the system by integration over all one- dimensional space. This 
involves, in particular, choosing a representation for the 7-matrices. Here we 
have to be careful, since a bad choice of representation can spoil manifest 
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reflection invariance (see Appendix A of Ref.|]I|). We choose the Pauli-Dirac 
representation, namely 



7o = 0-3 ; 7i = i°i and 75 = 7071 = o x . 



(7) 



Substituting the CGNM Hamiltonian written in Nambu basis in the dy- 
namical equation (|6|), we obtain an explicit dynamical equation which de- 
scribes the time evolution of the Nambu parameters. The calculation of 
traces is lengthy but straightforward. Taking the case N = 1 for simplicity 
and splitting the resultant complex equation into real and imaginary parts 
we have 
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where ii and I2 are the divergent integrals 
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A renormalization procedure is now required to deal with these diver- 
gences. In general, renormalization procedures consist in combining diver- 
gent terms with the bare mass and coupling constants of the theory to define 
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new, finite (or renormalized) values of these quantities. In other words, the 
bare mass and coupling constants are chosen to be cut-off dependent in a way 
that will cancel the divergences. In the present case, however the divergent 
integrals @ involve the dynamical variables themselves in the integrand, so 
that even their degree of divergence is not directly computable. In order to 
handle this situation we will use a self-consistent renormalization procedure 
inspired in Ref.||. 

The renormalization prescription that we use is based on the consideration 
of the static solutions of the dynamical equations (§), which satisfy 



We will show explicitly that controlling the divergences of these equations 
will also control the divergences which appear in the kinetic regime of the 
mean-field approximation [see below, Eqs.(|T7|)]. 

In order to obtain the renormalization prescription we introduce a regu- 
larizing momentum cut-off A and begin by assuming that, in order to render 
the theory finite, the bare coupling constant g 2 must approach zero for large 
values of the A as (see e.g. Refs.0, |J) 
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where the form of the first factor is dictated by later convenience. We next 
assume that the integrals I\ and I2 have logarithmic divergences 
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where a, b, c and d are finite constants. Substituting (|T|) and the ansatze 
fll3l) in static equations (III]) and (|TT1), we obtain 



d sin 7 k 
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tan2v9 k | eq = 
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We now verify that the assumed divergent character of the integrals I\ and 1% 
is consistent with the ansatze (0). Substituting the solution flU]) into Eq.(^) 
we find that Ji and I2 indeed have the prescribed logarithmic divergence (see 
Appendix C of Ref. Q). Moreover, from this calculation, we obtain the values 
of the constants a, b, c and d. We find 6 = 1, while d remains arbitrary. The 
renormalized static solution of our system in the mean-field approximation 
are then obtained by simply substituting these values in Eq.(|H|): 



dsin7 k | eq = 

(15) 

(-l) n m|k|rf 
tan2^ k | eq=[k2 + (i + d)m2] . 

We also show in Ref.jl]] that the connection between particle mass m and 
quasi-particle mass m eff is given by 

m e ff = (1 + d)m . (16) 

From the redefinition of the mass scale given by Eq . (|16D , we note that, un- 
like the situation found in connection with the 1/N expansion, the use of 
the Gaussian ansatz, Eq.(f|), parametrized by the canonical transformation 
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leading to the quasi-fermion basis, allows for the direct dynamical determina- 
tion of the stable equilibrium situation of the system [see Eq.(^)], including 
symmetry breaking [when d 7^ — 1, see also Ref. JTJ] ] and mass generation [see 
Eq.fllEDI- Moreover, the renormalization procedure effectively replaces the 
dimensionless coupling constant g 2 by the free parameter d associated to the 
mass scale [see Eq.flT6|)l. This is analogous to the phenomenon of dimensional 
transmutation found by Gross and Neveu || in the case of an 1/N expansion. 
Finally, aside from the over-all mass scale (characterized by d) there are no 
free adjustable parameters. 

Using Eqs.(||), (|T2| ) and fll3|), we finally write the renormalized form of 
the dynamical equations that describe the mean-field time evolution of the 
system. As mentioned before they are now also finite and read 

z>k,i = e z> k ,2 = 

/ x ,l k l 
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2.1 Small oscillations regime 

In order to study the small oscillation regime of the kinetic equations we 
next linearize the above kinetic equations around the static solution (|15|) 
taking = = 0. We begin by introducing the displacement away from 
equilibrium of the dynamical variables <£> k and 7k 

Vk = V^kleq + 5^k 

(18) 

7k = 7k|eq + #7k , 
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where the static solution y3 k | eq an d 7k|eq are obtained from the Eq.flTo]) [see 
also Ref. ill 



m|k|<i 
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The quantities Sip^ and 5^ ,will be treated as (first-order) small displace- 
ments. Functions of the dynamical variables are expanded also to first-order 
around the equilibrium solution (|19D. Therefore, we must linearize the diver- 
gent integrals (Qj around equilibrium. Taking z/ kA = we have 

h = /{ 0) + / 1 (1) + o(^ k ) 2 

dk' „ , „ rdk' 
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The linearized form of kinetic equations for Sip-^ and #7 k are then obtained 

as 



50k = r^Tk (21) 
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(22) 

where the renormalization procedure (|T2"|) controls the logarithmic divergence 
of the integral appearing in (|^) [see below Eqs. (|38|) and (^2|)1- Substituting 
Eq.(|22|) into Eq.(pT|) we obtain finally 

<5</3 k + 4[k 2 + (1 + d) 2 m 2 }5ip k + 

(23) 

-4 g) (* + DM / ^ [(k/)2+( j k ;', )2m2]1/2 ^ ^ 0. 



As usual in small oscillation treatments, this is a linear oscillator equation. 
Note that the last term couples different momenta. The solution to this 
problem involves determining the normal modes of small oscillation and their 
frequencies. This is done by looking for solutions of Eqs.(^Il) and (|22]) , which 
are of the form 



(24) 

57k = Ike*"* , 



where \I/ k and T k are time-independent amplitudes. Substituting the solution 
(33) into Eqs.(^Tp and (^), we obtain equations for these amplitudes, namely 



iujty k + — = (25) 
ko 

zo;^p^r k - 4[k 2 + (1 + dfm 2 }^^ + 

(26) 
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+4 g) ({ + l)|k| /rfk- [(k , )2 + ( J^ )2m2]1/2 ^ = 



Two-body dynamics from the linearized mean- 
field equations 



In this section we implement and discuss a reinterpretation of the small 
oscillations problem along the line proposed by Kerman and Lin in Ref.O. 
We begin rewriting equations (|25D and (|6|) as 



(kf) 2 ( g 2 \ , N Ikl r „ , Ik'l cj 2 



^--V^ (e+1) #)i^^ k ' = ^ k ' (28) 



where (kf) 2 = [k 2 + (1 + d) 2 m 2 } = [k 2 + m 2 ff ]. 

The crucial point is to realize that Eq. has the form of a Lippmann- 
Schwinger equation with separable potential 



(k\V\k') = v(k)v*(k') = (f) (£ + l)h(k)h(k') 
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In order to interpret this scattering problem we follow Ref.|| in relating 
dynamical small amplitude distortions of the Gaussian vacuum to two-quasi- 
fermion states. The Gaussian vacuum |0) can be written explicitly in terms 
of the fermion operators fej^ and the Nambu parameters in the well-known 
Bardeen-Cooper-Schriffer (BCS) form 
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|0> = I] [cosy? k + (-l^sin^e-^&J^ 

k>0,A^A 
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where y? k = y2 k | eq and 7 k = 7 k | eq and |0) is the vacuum of the 5-operators. 
It is also well known that |0) is the quasi-fermion vacuum, namely 



Ac,a|o) = o 

and that one- and two-quasi-fermion states can be written as 
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On the other hand, making the variations Eq. (]24) in Eq.(j30l) leads to a 
first order variation of the Gaussian vacuum which, appart from the expo- 
nential time dependence, is given by 



*|6> = vI/ kl 4 iiA ^l i)Ai |0)-^(-l) Ai r kl sin¥, kl | cq 6 kiiAi 6 kiiAi 



t 
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(34) 



This result just illustrates the well-known theorem by Thouless ]14[ when 
Tki = (as is appropriate in a static context), in which case it says that 
the linear variation corresponds to two-paired quasi-fermion addition to the 
Gaussian vacuum. Kerman and Lin have used this fact, in the context of 
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the 4 theory, to associate the scattering problem (|28| ) to the two-boson 
dynamics. A similar association can thus also be made here. Moreover, the 
last term of Eq. fl3"4] ) indicates the presence of dynamical symmetry breaking 
effects, since this term, proportional to the variation of the phase 7k, vanishes 
unless the chiral symmetry is broken in the static vacuum 7^ 0). 

The small oscilation regime can thus be seen as a non-perturbative way 
of approaching the dynamics of paired two-quasi-fermion excitations of the 
vacuum, including dynamical symmetry breaking effects. 

3.1 Analytical solution of the linearized equations 



We will now show how Eq.(|27|) and pq ) can be solved analytically. A general 
solution to two-fermion wave function will have two terms. The first one 
is the free solution (g = vanishing potential) and represents an incident 
wave. The second term is the non trivial part (when g ^ 0) which couples 
different momenta, and is associated with the scattered wave. Thus 



|k| 

Lett 
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tf(k,q;u;) = a<5(q - k) + 




(35) 

k'| 

7 *(k',q;u;) , 



where q is interpreted as the relative momentum for two incident quasi- 
fermions and a is an overall phase factor. We choose the outgoing wave 
condition (+ie) as solution of Eq. (p8[) , but we could have chosen e.g. the 



incoming wave condition (— ie) or Van Kampen wave condition |15| or another 
condition. 

Integrating the Eq.(|35|) with respect to k 



dk TjA*^* w = 7 — 2 T (36) 

( k o ) h - (o2/4tt) (£ + 1) _2e^ 



and substituting this result back into Eq.(]3q) yields a general solution for 
*(k,q;w) 
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where A + (a;) is given by 
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The oscillation amplitude r(k, q; cu) is obtained from the Eqs. (p7|) and 
and reads 



(39) 
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kf (\k\\ 1 /|k| 
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Finally, substituting Eq.(]3"7D into Eq.fl2"E|) we obtain the oscillation fre- 
quencies 

u = 2qf = 2[q 2 + m 2 ^} 1 / 2 , (40) 

where q is the relative momentum for two incident quasi-fermion with masses 
m e ff = (1 + d)m. 

We observe that we can understand the factor 2 in the frequencies of 
oscillation u [see Eq.([H])], as related to the treatment of harmonic oscilla- 
tors in terms of the sympletic groups given by Goshen and Lipkin It 



can be interpreted classically by noticing that, since for harmonic oscillators 
the frequency does not depend on the amplitude of the motion, if a set of 
independent particles in a harmonic field is symmetrically stretched out of 
equilibrium, it will subsequently pulsate with frequency 2u, where u is the 
frequency of oscillation of the independent particles. 
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4 Bound states from the small oscillations 
regime 



In this section we will examine the condition for existence of bound states in 
the small oscillation regime around the stationary solution (vacuum) of our 
fermionic system. 

From Eqs.fl28|) and (p9|) we verify that the potential term, which describes 
the time evolution of our system in this regime is separable. Again, in analogy 
with scattering theory, we can evaluate the corresponding T matrix JF]. We 
find 

T(k,k>)oc/i(k)^-^(k') (41) 



with h(k) given by Eq.(|2T3|) and A + (u) given by Eq.(^). The bound states 
are given by the poles of the T matrix. Therefore, we search for the zeros 
of A + (u). It is clear that the integral in A + (u) contains a logarithmic di- 
vergence. To keep it under control, we use the renormalization procedure of 
the coupling constant given by Eq.([[|). Substituting Eq . (|i~2"D into Eq.(|38|) 
we get 



{UJ) ~ 11 [m* J J -a [k 2 + m c 2 ff ]V2 [k 2 + m 2 ff - oo 2 /A + ie] ' ( } 



where we introduce the regularizing momentum cut-off A. 

In the interval < uj < 2m e g the integral of A + (u;) is well defined and 
we can set e = 0. A straightforward calculation yields 

A+(w) = 2[/(w)-j(d)] (43) 



where f(u>) 



"4m 2 ff " 


1/2 




"4m c 2 ff " 


-1/2 ^ 


arctan < 












. ^ 2 





16 



(44) 

m = ln [\ITd\_ ■ 

Fig.l shows the zero of the A + (u) as function of d. In this calculation 
q = 0, therefore to is the mass of the bound state. Obviously, when (1+d) = 
or m cff = (free system, see Ref.fi]) there is no bound state. We see from 
Fig.l that a bound state of quasi-fermions occurs when 0.74 < (1 + d) < 2, 
and that the mass of this bound state will vary in the interval < u < 2m e g. 
Gross and Neveu obtain M a = 2M F for the mass of the a particle in 
leading- 1/N aproximation, where Mp is equivalent to m e ff. They argue that 
in higher order they might find that 

M (T = 2M F [1 + 0(1/N)} . 

From Fig.l we verify that u = 2m e g corresponds to (1 + d) =2. Observing 
that j(l+d = 2) = 0, we may conclude that j(d) can be contribuition 
of higher order to the Gross- Neveu result. 

We believe that in the limit N — > oo the function j(d) — > 0. On the other 
hand, when iV is finite [N = 1 in our calculation - see Eq.®], the mass u of 
the bound state depends of the renormalized coupling constant d as shown 
in Fig.l. This dependence can not be obtained from 1/N approximation. 

Therefore, we can conclure that to N —>■ oo the 1/N approximation and 
our mean-field approximation are equivalent. On the other hand, when iV is 
finite our approximation permits to obtain the higher order contribution to 
the Gross- Neveu result [§]. 

It is important to observe that the higher order term obtained in Eqs. ( |43| ) 
and (0) in this approach does not contain all terms of order 1/N, since the 
mean-field approximation is not a systematic expansion in the parameter 
1/JV. 

Surprisingly, we have obtained for the function A + (u) a structure which 
entirely reproduces that which has been found by Kerman and Lin in 
their study of the bosonic A0 4 theory in terms of a Gaussian time-dependent 
variational approach. 

Finally, when uj > 2m c g the integrand of A + (uj) has a singularity at 

k = ±^Juj 2 /A — . From the theory of residues we obtain 
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l + (l-4m e Va; 2 )^ 
l-(l-4my^) 1/2 J 

(45) 
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Now A + (c<j) does not have any zeros. The interesting point here is to 
observe that 

lim A + (uj) — ► In 

u>— >oo 

so that in the large frequency limit (u — > oo) the T matrix goes asymptot- 
ically to zero. We thus recover, in the present approximation, the assymp- 
totically free character of the CGNM. 



™eff, 



OO 



(46) 



5 Discussion and conclusions 



In Ref.ffl we described a treatment of the initial- values problem in a quan- 
tum field theory of self-interacting fermions in the Gaussian approximation. 
Although the procedure is quite general, we implemented it for the vacuum of 
an uniform (1+1) dimensional relativistic many-fermion system described by 
Chiral Gross-Neveu model (CGNM). We obtained the renormalized kinetic 
equations which describe the effective dynamics of the Gaussian observables 
in the mean-field approximation for this system. 

In this work, we have considered the linearized form of the mean-field 
kinetic equations obtained in Ref. U around the stationary (vacuum) solu- 
tion. The two-quasi-fermion physics can be analytically investigated in this 
approach. In particular, we have solved these equations completely. From 
the solutions, we have reinterpreted the near equilibrium physics of our sys- 
tem as a problem of quasi-fermion scattering and have found the condition 
for the existence of a quasi-fermion bound state. 
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We verify that for N finite (in this work N = 1), the bound state mass 
obtained from our approach contains a term which depends of the renormal- 
ized coupling constant as can be seen in the Fig.l. In the case of an 1/N 
expansion p] this dependence can not be found, so in the limit N — > oo this 
term goes to zero. Therefore, to small N, our approach permits to obtain 
the higher order contribuition to the 1/N expansion. 

Finally, it is important to observe that the higher order term obtained to 
the bound state mass from our approach contains no necessarily all terms of 
1/N order, since the mean-field approximation is not an 1/N expansion. 
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Figure Captions 

Figure 1 - The curve represents the mass u of the two-quasi-fermion 
bound state in small oscillation regime for our system as a function of the 
renormalized coupling constant d. 
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This figure "figl-l.png" is available in "png" format from: 
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